Anharmonic suppression of Charge density wave in 2H-NbS 2 
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The temperature dependence of the phonon spectrum in the superconducting transition metal 
dichalcogenide 2H-NbS2 is measured by diffuse and inelastic x-ray scattering. A deep, wide and 
strongly temperature dependent softening, of the two lowest energy longitudinal phonons bands, 
appears along the TM symmetry line in reciprocal space. In sharp contrast to the iso-electronic 
compounds 2H-NbSe2, the soft phonons energies are finite, even at very low temperature, and no 
charge density wave instability occurs, in disagreement with harmonic ab-initio calculations. We 
show that 2H-NbS2 is at the verge of the charge density wave transition and its occurrence is only 
suppressed by the large anharmonic effects. Moreover, the anharmonicity and the electron phonon 
coupling both show a strong in-plane anisotropy. 

PACS numbers: 78.70.Ck,74.25.Kc,63.20.Ry,71.45.Lr 



I. INTRODUCTION 

In transition metal dichalcogenides (TMD) the in- 
terplay between strong electron-electron and electron- 
phonon interactions produces a wide variety of instabil- 
ities ranging from charge density wave (CDW), as in 
TaSe 2 , TaS 2 , NbSe 2 and TiSeJ™, to Mott-insulating 
states in TaS 2 4 . The most eminent representative of this 
class of materials is 2H-NbSe 2 in which a CDW (Xcdw = 
33 K) coexists with superconductivity (Tq = 7.1 K) with 
a strongly anisotropic superconducting gap. The debate 
about the origin of the CDW instability in this system 
has recently been revived. Several independent studies 
tried to locate the CDW gap using Angle Resolved Photo 
Emission Spectroscopy (ARPESp^. These studies how- 
ever found either Fermi surface nesting inconsistent with 
the CDW wave vector or no Fermi surface nesting at 
all. As such these studies were not conclusive concern- 
ing the origin of the charge density wave instability in 
TMD. On the other hand, theoretical calculation ^ Q * n * 
based on density functional theory (DFT) reproduce the 
experimental CDW state, despite a real part of the bare 
electronic susceptibility having only a broad maximum 
rather than a sharp divergence. In contradiction with 
a pure nesting scenario, the authors conclude that the 
main driving force for the CDW transition is the electron- 
phonon coupling (EPC). The accuracy of these theoret- 
ical calculations was recently validated against experi- 
ment by direct measurements of the phonon dispersion 
using inelastic x-ray scattering (IXS/^. 

In most metallic TMD, superconductivity coexists 



with a charge ordered state or at least, it occurs in prox- 
imity of such a state. For many years it has been thought 
that the occurrence of a CDW state was mandatory for 
superconductivity. However there exists at least one case 
that disagrees with this picture, namely 2H-NbS 2 . This 
compound is isoelectronic and isostructural to 2H-NbSe 2 , 
and has a similar superconducting transition tempera- 
ture ( Tq = 6 K) . However no CDW order occurs in 2H- 
NbS 2 ^^, which does not fit really well with the variety 
of electronic instabilities found in the TMD family. 



In this work we measure the phonon dispersion and 
linewidth of 2H-NbS 2 , using thermal diffuse scattering 
(TDS) and IXS. Moreover, we perform first-principles 
linear-response calculations of the harmonic phonon dis- 
persion. Surprisingly we find that while theoretical cal- 
culations do predict the occurrence of a CDW state, in 
experiments only a wide and marked softening of the 
low energy phonon branches is present, meaning that no 
CDW state occurs down to 2K. We demonstrate that 
this disagreement is due to the occurrence of large anhar- 
monic effects that prevent CDW formation and stabilize 
the lattice. Our work shows that 2H-NbS 2 is, actually, al- 
most as close to the CDW instability as all other metallic 
TMD and points out the relevance of anharmonic effects 
to describe the competition between CDW and supercon- 
ductivity in TMD. Moreover, we show that the EPC is 
strongly anisotropic, giving a natural explanation to the 
superconducting gap anisotropy. 
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FIG. 1. Reconstruction of the x-ray diffraction in the (H,K,0) 
plane of 2H-NbS 2 at 300 K and 77 K, from a 3D mapping 
of thermal diffuse scattering. A symmetrization has been ap- 
plied. Diffuse intensity is present only along TM and, more 
precisely, only where scattering geometry selects longitudinal 
phonons (white ellipses). No diffuse features are visible for 
transverse geometry (black ellipse) . At 77K small extra spots 
are visible. 



II. EXPERIMENTAL DETAILS 

High quality single crystals of 2H-NbS2 were prepared 
by iodine vapor transport with a large excess of sulphur 
in a sealed quartz tube. Powder diffraction on several 
batches gave a predominance of 99% for 2H polytype 
(P6s/mmc space group), and 1% for 3R. The polytype 
was then fully checked by x-ray diffraction for each se- 
lected single crystal. Refined lattice parameters at room 
temperature are a = b = 3.33 A, c = 11.95 A. The super- 
conducting properties of the main s ample used for this 
study were published elsewher e 1 15 * 16 *. 

TDS experiments 17 at 77 K and room temperature 
were performed using a PILATUS 6M detector and a 
wavelength A = 0.6966 A (Fig.[l]). In the (a,b) plane, the 
Bragg peaks are sharp: full width of the (1,1,0) reflection 
is 0.13°, evidencing a good structural order. Meanwhile 
along the c-axis, analysis of the diffraction pattern re- 
vealed a correct 2H stacking over 3 unit cells on average, 
in agreement with a previous study^. Between the Bragg 
spots in the (a,b) plane, diffuse lines can be detected at 
room temperature, and their intensities increase at low 
temperature. These lines are present only along the TM 
direction and, more precisely, only where scattering ge- 
ometry selects longitudinal phonons (circled in white in 
Fig. [I]). No such diffuse features are visible for transverse 
polarization (circled in black in Fig. [T]). This indicates 
that an anomaly occurs for longitudinal phonons along 
TM. In addition, at low temperature, we observe the ap- 



parition of tiny spots at wavevectors q = (0.38,0.16,0) 
and at the M point q = (0.5, 0, 0). 

To get further insight regarding these diffuse features, 
we performed IXS measurements with a photon energy 
of 17.794 keV using the (9,9,9) reflection on the high- 
resolution silicon backscattering monochromator. Cor- 
responding instrumental energy and momentum resolu- 
tions were 2.6 meV FWHM (determined from the least- 
square fit of the Lorentzian curve of the elastic peak), 
0.017 A" 1 in the (H,K,0) plane and 0.04 A" 1 in the 
(0,0, L) direction, respectively. The x-ray beam was fo- 
cused down to 200 x 60 fim (width x height) on a 
single crystal of 2H-NbS 2 of 450 2 x 100 fim (a x b x 
c). We found the tiny spots to be very weak elastic 
peaks, considering the high electronic susceptibility at 
q = (0.5,0,0), these may be traces of a charge density 
modulation induced by local strain or of Friedel's oscil- 
lations. Meanwhile, the diffuse features were due to a 
strong phonon softening. This can be evidenced when 
mapping the phonon dispersion (see Fig. [2]) along the 
TK and TM line at fourteen different temperatures be- 
tween 300 K and 2 K (pumped 4 He cryostat). We per- 
formed our measurements close to the (3,0,0) Bragg re- 
flection, where we should observe both supposedly soft 
modes on LO and LA branches, according to structure 
factor calculations. The presence of these two modes 
is confirmed experimentally by the asymmetric profile 
of the experimental spectra in Fig. [3] (b). Harmonic 
phonon frequencies were calculated in the framework of 
DFT in the linear response^. We use the local den- 
sity approximation. The dynamical matrix were cal- 
culated on a 10 x 10 x 1 phonon-momentum grid in 
the Brillouin zone and Fourier interpolated throughout 
the Brillouin zone. The electronic integration was per- 
formed using 20 x 20 x 6 electron-momentum grid and 
an Hermite- Gaussian smearing (electronic temperature) 
T e i ec = 0.136 eV. The EPC contribution to the phonon 
linewidth (a quantity that does not depend on real or 
imaginary nature of the phonon frequency) has also 
been calculated following RefP4 Experimental spectra 
were fitted using standard 21 damped harmonic oscillator 
(DHO) functions for the phonons with the Bose factor 
included. The linewidths were constrained to calculated 
electron-phonon ones (see Tab. |l|, convoluted by the 
quasi-Lorentzian experimental resolution function. This 
procedure yields excellent results as seen in Fig. [3j and 
the resulting experimental dispersion of the phonons is 
plotted in Fig. [2j together with the calculated harmonic 
phonon frequencies. 



III. RESULTS AND DISCUSSION 

Along the TK direction the low energy phonon dis- 
persion is in excellent agreement with harmonic first- 
principles calculations, as also shown by the scans in 
Fig. [3] (a). On the contrary along TM and close to the M 
point, harmonic calculations do show imaginary phonon 
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FIG. 2. (Color online) Experimental phonon spectrum of 2H- 
NbS2 at room temperature (full red circles) and 2K (open blue 
diamond), compared to ab initio, zero temperature, harmonic 
phonons calculation (dashed lines). The two soft phonon 
branches along TM are nearly degenerate at room temper- 
ature, and slightly split away as temperature decreases. At 
the exception of the scans along TM at room temperature, we 
limited ourselves to energies below 30 meV. 



frequencies for the lowest energy phonon modes, mean- 
ing that the high temperature structure becomes unsta- 
ble and that a CDW with ordering vector close to the 
0.5 a* occurs. Experiments find a completely different 
picture. Namely, a marked softening of these branches 
occurs but the system remains dynamically stable with 
no indication of CDW instabilities. 

More insight can be obtained by analyzing the be- 
haviour of the measured phonon frequencies as a function 
of temperature. Experimentally, along the TK direction 
(see Fig. |2| the phonon frequencies are essentially tem- 
perature independent. On the contrary, along TM, and 
in particular close to the M point, the low energy phonon 
modes are strongly temperature dependent. This is bet- 
ter seen in Fig.[4]where the lowest energy phonon frequen- 
cies at 0.75 TM are plotted as a function of temperature. 
The phonon frequencies are more than doubled in the 
2 — 300 K temperature range, a behaviour that can only 
be explained by invoking strong anharmonic effects^, as 
discussed hereafter. 

In solids, anharmonicity has two main effects. Firstly, 
it is the origin of the thermal expansion in solids. How- 
ever this contribution is negligible here, as the lattice 
parameters barely change from zero to room tempera- 
ture (a(300K) = 3.3295 A, a(2K) = 3.3230 A). Secondly, 
when the lattice parameters are fixed, phonon-phonon 
scattering results in a temperature dependent phonon- 
frequency shift to higher energies and in a tempera- 
ture dependent enhancement of the phonon-linewidtfP^. 
These contributions add to the normal electron-phonon 
contributions that are, on the contrary, temperature in- 
dependent. Thus, a strong temperature dependence of 
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FIG. 3. (Color online) a) IXS scans at (H,K,L) = (3.184,- 
0.072,0) fitted by the convolution of a DHO (with Bose factor) 
and the Lorentzian experimental resolution. Ab initio phonon 
calculations, convoluted by the experimental resolution, are 
in good agreement with experimental results, 
b) IXS scans at (H,K,L) = (3.375,0,0) for different tempera- 
tures. For clarity, scans are shifted up. In the superimposed 
fits, the linewidths and amplitude ratio of the two phonons 
are fixed according to the ab-initio calculations (see text). 
Ab initio phonon calculations, convoluted by the experimen- 
tal resolution, are also in good agreement with experimental 
results when using higher T e i ec - 



the phonon frequencies is a fingerprint of large anhar- 
monic effects. In fact, we find here that the square of 
the two phonons frequencies at the softest q vector in- 
creases linearly with temperature (Fig. [4|. This power 
law dependence of the frequency, also observed in 2H- 
NbSe2 12 , is typical for a soft mode. As temperature is 
increased, the unstable harmonic oscillator is stabilized 
by higher order anharmonic potential, induced by the 
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FIG. 4. Experimental phonon energies, squared, from our 
fitting procedure with two phonons (see text), as a function of 
temperature at (h,k,l)=(3.375,0,0). Energies are compatible 
with the mean- field theory which predicts a power law T2. 
Lines are guides to the eyes. 

phonon-phonon interactional The critical exponent 1/2 
indicates that the phonon-phonon interaction can be de- 
scribed in a mean- field approach, i. e. considering the soft 
mode in equilibrium with the average field induced by all 
other phonons 2 ^. 

First principles evaluation of third and fourth or- 
der perturbative coefficients in phonon-phonon scattering 
cannot be easily done here as the harmonic solution has 
imaginary phonon frequencies as seen in Fig. [2] As such, 
the anharmonic correction is as big as the bare harmonic 
phonon frequency itself. Non perturbative approaches 
which include anharmonic effects 2 -^ are required. These 
approaches are however very time consuming and unfea- 
sible for the case of 2H-NbS2 with 6 atoms per cell. In 
order to interpret the IXS spectra along TM, we artifi- 
cially smear out the sharp Fermi surface with an elec- 
tronic temperature a technique also used in RefP. 

We have calculated the contribution of the electron- 
phonon coupling to the linewidths of the phonon modes, 
along the two high symmetry lines YM and TK. The val- 
ues of the half width at half maximum are shown in ta- 
ble [T] Along the TK direction the linewidths of the first six 
branches are far below the experimental resolution (1.3 
meV). However, along the TM direction, the electron- 
phonon coupling is responsible of a large broadening of 
the two soft modes on a wide range of wavevectors. In 
both directions, the linewidth is almost independent of 
the smearing of the Fermi surface, driven artificially by 
the increase of T e ^ ec , up to T e i ec = 0.34 eV and consis- 
tently with previous reportP^ (Fig.4(b)). 

On the contrary, T e i ec has a very strong effect on the 
renormalisation of the phonon energies close to the M 
point. This is illustrated in Fig. [5] Yet, a negligible 
renormalization is observed along the TK direction, thus 
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TABLE I. Phonons linewidths 7 in meV (Half Width Half 
Maximum) of the first 6 phonon modes, calculated ab initio 
for several positions along TM (top) and TK (bottom). Posi- 
tions are expressed as a fraction x of TM (resp.TK). At each 
position x, the 7 are sorted by increasing order of the phonon 
energy. 



emphasizing the anisotropy of the electron-phonon cou- 
pling. 

Typical IXS spectra calculated at 0.75TM for elec- 
tronic temperatures of T e / ec =0.34 eV and 0.27 eV (in- 
cluding the electron-phonon linewidth, weakly dependent 
on the chosen electronic temperature, as well as the ex- 
perimental resolution) are shown in Fig. [3] (b) and are in 
good agreement with the experimental data. The elec- 
tronic temperature smears the Fermi surface, and by do- 
ing so suppresses the effect of the EPC. This in turn 
increases the phonon frequency artificially and stabilizes 
the lattice. However, this increase of T e \ ec is not the cor- 
rect physical picture since to reproduce changes in spec- 
tra in the 0-300 K temperature range, the electronic tem- 
perature must be changed by at least 0.1 eV. and even 
more in the case of 2H-NbSe^ or lT-TiSej^. Thus, 
considering this effect as the origin of the tem peratu re 
behaviour of the phonon frequencies is incorrect l 12 l 25 l. 

On the contrary, the aforementioned anharmonicity, 
neglected in these first-principle calculations, is likely 
responsible for these discrepancies. This was already 
suggested by Varma and Simons 2 ^, who demonstrated 
that anharmonicity can be responsible for a reduction 
by an order of magnitude of Tqdw m lT-TaS2, hence 
accounting for the anomalously large ratio of the CDW 
gap amplitude to Tqdw m these compounds. The strong 
mean field power law temperature dependence of the soft 
phonon modes of 2H-NbS2, and the prediction by har- 
monic DFT calculations of a CDW instability close to 
the M point, both consistently indicate that anharmonic 
effects are, in this case, so strong that they actually sup- 
press the CDW and stabilize the crystal structure, leav- 
ing however the superconducting properties essentially 
unaffected 2 ^. 

In real space, the main components of the soft phonon 
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FIG. 5. Calculated phonon energy squared for the two lowest 
branches at 0.8rM as a function of the smearing of the Fermi 
surface with T e i ec - 



modes are the vibrations of the Nb atoms in the (a,b) 
plane, but also the vibrations along the c-axis of a cou- 
ple of S atoms situated just above and below the plane of 
Nb atoms. The longitudinal displacement of Nb atoms 
along TM corresponds to a movement toward this couple 
of S atoms. Interestingly, the distance between these S 
atoms is unusually short (2.97 A,while typical S-S bond 
are 3.3-3.4 A^l). It therefore suggests that the displace- 
ment of Nb atoms is hampered by the great elastic cost 
of moving apart the S atoms, which would be the micro- 
scopic origin of the anharmonicity. The decisive role of 
the force constant between the S atoms in the same S-Nb- 
S sandwich was already emphasized in a study based on 
a tight binding approach 29 and with first principles cal- 
culations in TiSe2 30 . In this respect 2H-NbS2 is rather 
unique as it is on the verge of the CDW transition and 
is only stabilized by anharmonic effects. 

A second important result of our work is that the EPC 
(responsible for the softening close to the M point) is 
strongly anisotropic throughout the Brillouin zone. The 
ab-initio calculation yields a very large phonon linewidth 
of 1.5 meV at 0.8 TM, due to the EPC, compared to less 
than 0.1 meV along the TK line. This strong in-plane 
anisotropy has necessarily important consequences for 
the superconducting properties of TMD. In this respect it 
is instructive to compare superconducting gap measure- 



ments of the two isoelectronic dichalcogenides 2H-NbSe2 
and 2H-NbS2- In the former, there exists two main val- 
ues of the superconducting gap EHHUl This was initially 
attributed to the CDW order, pre sumably competing 
with the superconducting instabilit ^ 7 | 9 | 36 | 37 l However a 
similar observation in the superconducting state of 2H- 
NbS2 by the same various methods (scanning tunneling 
spectroscopy, specific heat measurements, H&, magnetic 
penetratio n depth m easurement), clearly ruled out this 
assumptior l 13 * 15 * 1 ^^. As a consequence of the strongly 
anisotropic softening of the phonons that we report in 
this study, the EPC, and therefore the pairing strength, 
are themselves highly anisotropic. This suggests that 
the two superconducting gap amplitudes observed in 2H- 
NbS2, and by extension in 2H-NbSe2, rather originate 
from the anisotropy of the electron-phonon coupling. 



IV. CONCLUSION 

In conclusion, we identified two wide soft phonon 
modes in 2H-NbS2- They strongly soften along the TM 
direction as the temperature is lowered, but their energies 
always remain finite and never freeze into a static mod- 
ulation of the lattice. This absence of a CDW is in dis- 
agreement with the first-principle harmonic calculations. 
The mean field power law temperature dependence over 
the range 300 K-2 K is also an unambiguous evidence of 
strong anharmonic effects which smear the effect of the 
EPC. These anharmonic effects are responsible for the 
absence of CDW in 2H-NbS2, but also for the low Tqdw 
transition temperatures in other TMD. Thus, 2H-NbS2, 
which is the lightest compound among the TMD, is a 
TMD on the verge of a CDW, only hampered by anhar- 
monic effects. These, however, leave the superconducting 
properties essentially unaffected compared to 2H-NbSe2, 
which suggests that the relevant parameters for super- 
conductivity in TMD are the anisotropy and strength of 
the electron-phonon coupling, rather than the nature of 
the ground state. 
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